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Abstract. The critical behaviour of directed self-avoiding walks is studied on 
parabolic-like systems with a free boundary at x = ±Ct a . Using a scaling argument, 
1/C7 is shown to be a marginal variable when a = v_i/ u \\ =1/2 i.e. on a parabola. 
As a consequence the directed walk may display varying local exponents. Such a 
behaviour is indeed observed for restricted walks. This generalizes a result of Cardy 
showing that nonuniversal behaviour occurs at corners for isotropic systems. 



Isotropic systems are known to display nonuniversal critical behaviour at corners 
[1-6], local exponents varying continuously with the opening angle 9. As shown 
recently [7], this may be linked to the scale invariance of these shapes in the case 
of isotropic systems. In parabolic-like geometries where the boundary is located 
at X(t) = ±Ct a an isotropic change of scale transforms C into b OL ~ 1 C where b is 
the dilatation factor. The dimension of C vanishes when a = 1, i.e. in the corner 
geometry, then C = tan#/2 is a marginal variable leading to ^-dependent exponents. 
When a > 1, C grows under renormalization and the critical behaviour is that of a flat 
surface whereas when a<l, C decreases and one gets either a line geometry or a cut, 
depending on the location of the system relative to its boundary. 

These considerations may be extended to the case of anisotropic systems [8] 
for which the correlation length diverges at the critical point with different critical 
exponents in the parallel and perpendicular directions. If the t-axis of the boundary 
defined above coincides with the parallel direction of the system with a correlation 
length exponent z/y = Zfj_, under an anisotropic change of scale [9] with dilatation 
factors 6|| = b z and b ± =b, C is changed into 

C' = b za ~ x C (1) 

and marginal behaviour is obtained when a = l/z. 

In the present work, we check these ideas in the case of the directed self-avoiding 
walk [8,10] for which z = 2 so that varying exponents are expected inside a parabola. 
Due to the directedness, one cannot get any boundary effect for a walk outside a 
parabola. Although in the following we restrict ourselves to the 2D problem, similar 
results are expected in higher dimensions. 
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Let us consider a directed walk on a rectangular lattice in the (x, t) plane. At each 
step with At = r in the time direction, the walker performs a jump Ax = ±a towards 
one of the two nearest sites with the same probability so that the walk is restricted and 
may be also considered as a ID diffusion process. Furthemore the walker is assumed 
to start at the origin (x = t = 0) and remains confined inside a parabolic-like domain 
which for convenience we consider to be slightly shifted backwards in time 

\x(t)\ <X(t) =C(t + V ) a (2) 

The number of iV-step walks starting from the origin and reaching x = no at time 
t = Nr may be written as 

A^(0, no) = 2 N P(na, Nr) (3) 

where the front factor on the right gives the total number of unconfined walks with N 
steps and P(x, t) gives the probability to reach x = na without crossing the frontier. 
This probability satisfies the recursion equation 

P(x, t + t) = X - [P(x + a,t) + P(x - a, £)] P(x, 0) = 6 xfi (4) 

with the boundary condition 

P(x,t) = \x\=X(t) (5) 
In the continuum limit (a<Sl, r^Cl, a 2 /r = l), one gets the diffusion equation 

dP 1 d 2 P 



dt 2 dx 2 



(6) 



together with the time-dependent absorbing boundary condition given in (5). 

In a free random walk x 2 (t)=t so that when a>l/2 the walker no longer sees the 
frontier when t^>t* = (7 2 /( 1 ~ 2q and the probability distribution tends asymptotically 
to the Gaussian form, namely 

This is the behaviour expected from equation (1) since C grows under renormalization 
when a > 1 /2 so that the boundary evolves towards a flat surface geometry. 

When a < 1/2 the boundary changes the asymptotic behaviour and in order to 
work with constant boundary conditions, it is convenient to introduce the new variable 

leading to the following equation for P{y, t) 

f)P 1 f) 2 P BP 



with absorbing boundary conditions at y = ±C. 
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In the appendix, scaling arguments are used to get the form of the probability 
distribution. Its leading behaviour at long time when a < 1/2 is found to take the 
form 

r,/ x 1 ( \ f TTX \ 

— exp -— y- — cos — — 10 

v ' ; Ct a l \ 8C 2 I -2a) \2Ct a ) v ' 

This asymptotic expression becomes exact in the limit of a strip geometry when 
a = 0. A similar stretched exponential behaviour was previously observed in isotropic 
parabolic systems with relevant boundary effects [7]. 

In the marginal case a = 1/2 the variables separate in (9) and the problem is 
exactly solvable. The diffusion equation becomes 

.dP ld 2 P ydP 

and looking for P as a product 4>{i)ip{y) enables us to write down 

(* + "))^ = -A 2 ^ (12a) 



dV + n z )«L. XH , = o (12b) 



dz 2 \2 J dz 

where we used the new variable z = —y 2 /2 in the eigenvalue equation (12b). 
One recognizes Kummer's equation [11] so that the eigenfunctions are confluent 
hypergeometric functions. Even solutions 1 F 1 (A 2 ,^;^) must be selected since the 
initial condition introduces a reflection symmetry with respect to the time axis. These 
have the following z-expansion 

t-, / , \ ^ a a ( a + 1) z 2 

^ (a '^ ) = 1+ r + 4riy^ + -- 



a(a + l)(a + 2)---(a + n-l) z n , 
6(6 + 1)(6 + 2) • • • {b + n - 1) n\ [ } 

In order to satisfy the boundary conditions, the eigenvalues have to belong to 
a discrete set which, according to (2) and (8), corresponds to zeros of the confluent 
hypergeometric function for z = —C 2 /2. They are solutions of the implicit equation 

M>tn ,h-Tl= (14) 



2' 2 f 

The first differential equation (12a) simply yields 

<P m (t)~(t + V )- x ™ (15) 
and putting these together, one may finally write 



2 



(16) 
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Figure 1. Inverse of the largest eigenvalue of equation (12b) as a function of the 
square of the marginal parameter C. 



where the coefficients A m {rj) have to be chosen to satisfy the initial condition. The 
asymptotic behaviour is governed by the lowest eigenvalue Aq which was studied 
numerically (figure 1). 

Some analytical results may be obtained only in limiting cases. Let us first consider 
the strong curvature limit e = C 2 /2 — > 0. Then, according to (13) and (14) A^ is 
0{e~ 1 ) and introducing u = C 2 \ 2 n = 0(1) in (14), one gets cos(v / 2~u) =0(e) so that 
the eigenvalue spectrum reads 



K n = \ m+- 



i \ 2 2 
1 \ 7l Z 



+ 0(1) 



m = 0,1,2- •• 



2) 2C 2 

The eigenfunctions are obtained in the same way leading to 



(17) 



1-^1 ( Km o' 



X 



2' 2(t + rj) 
The completness relation 



= cos 



1 

m H — 



TlX 



2) C^t + v 



+ 0(C 2 ) 



(18) 



— — — cos 



1 A nx 
m + - 



2 J Cyft. 



8(x) 



(19) 



may be used to deduce the coefficients A m (rf) satisfying the initial condition P(x, 0) = 
8(x) and the probability distribution behaves as 



cos 



1 

m + — 



7TX 



2 J C/t + v 



+ 0(C) 



(20) 



the asymptotic behaviour being governed by the first term in the eigenvalue expansion, 
namely 



_ vr 2 f nx \ 

P(x,t)~t-^ cos ^) 



t > 1 



(21) 
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This is just the result given by (10) when a — > 1/2 if t 1_2a is changed into t l ~ 2ot — 1 
in the exponential, which amounts to modify a constant prefactor, in order to get a 
meaningfull limit. 

Let us now consider the opposite limit. First, when C is infinite, the boundary 
condition is satisfied with Aq = 1/2 since [11] 

aQ-H-*" (22) 

which corresponds to the Gaussian distribution 

(23) 



P(x,t) = — = exp 



2 

X 



2(t + V ) 



satisfying the initial condition when the limit r\ — > is taken. This is the free diffusion 
result (7) as expected for a walker starting on a surface which is flat in this limit. 
When the curvature is weak (C^>1) one may use the z~ ^expansion of the confluent 
hypergeometric function 

^(a, 6; z) = ^e z z a - b [l + O^ 1 )] + (_ z )-[i + (24) 

to find out the leading correction to 1/2. After some algebra one gets 

A » = 5 + ^ exp Hr) |1+0(£)1 (*) 

where e is the correction term itself. 

Let us finally turn to the evaluation of the critical exponents. Besides the radius of 
gyration exponents z/y =1 and zv_|_ = l/2, in analogy with thermal critical phenomena, 
one defines a bulk susceptibility exponent 7 which enters into the asymptotic behaviour 
of the total number of iV-step directed walks starting from the origin on an unlimited 
lattice, M N ~ n N N 1-1 [8]. Since for the restricted walk, N = t when the time is 
measured in r units, by (3) this can be written as 

/+00 
P(x,t)dx = 2 N (26) 
-00 

so that \i = 2 and 7 = 1. Surface exponents may be defined in a system with a 
straight free surface along the time axis by considering the total number of iV-step 
walks starting near the surface at x = 5 

/•OO 

/ Af N (5,x)dx ~ ^iV 71-1 (27) 
Jo 

and the number of iV-step walks starting and ending near the surface 

AUM) ~ /^iV 711-1 (28) 

Using the Gaussian distribution (7) and the method of images [10], it is easy to check 
that 7-l = 1/2 and 7 n = — 1/2 in agreement with the scaling law 2^ 1 —^ 11 =^ + u ± [8]. 
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Now, in the parabolic geometry, one may define two new exponents by considering 
either the total number of iV-step walks starting near the tip 

Af N (0,x)dx~ ^N^- 1 (29) 

'-Ct<* 

or the number of iV-step walks starting near the tip and ending near the boundary 

Af N (0, Ct a -8)~ fi N N 101-1 (30) 

Using the analogy with thermal critical phenomena, it may be shown that when 
l/2>a>0 these exponents satisfy the scaling law 

7o - 7oi = 7i - 7n (31) 
since the boundary is asymptotically flat and parallel to the time axis. When a>l/2 
one evidently gets 7 = 7 = 1 and 7 01 remains undefined since at long time the walks 
cannot reach the surface. In the marginal case, a = 1/2, the leading contribution to 
P(x, t) in (16) gives 

which combines with (3), (16) and (29) to yield 

7o = \ ~ Ao (33) 

The second exponent is obtained in the same way through a first-order expansion of 
the confluent hypergeo metric function. Then, by (14) 

/ l2 l C 2 C5\ 1/2 

and finally 

7oi = \ ~ A o (35) 

follows in agreement with the scaling law (31). When a < 1/2, by (3) and (10), the 
connective constant [i is given by 

ln// = ln2- lim — -- — (36) 

p n^oo 8C 2 l-2a v ' 

so that fx = 2 when a > and due to the exponential decay of the probability 

distribution, 7 = 7 01 = — oo. It follows that when C varies from zero to infinity, 

the marginal exponents interpolate between their values below and above a = 1/2. 

When a = 0, by (36), the connective constant is changed into 

^ = 2eXP ("^) (37) 

and as a consequence of the one-dimensionality of the system 7 = 7 01 = 1. Let us 
mention that doing the calculation with the transfer matrix technique on a lattice 
would give \i = 2 cos(7r/2(C + l)) instead of (37) but the continuum limit used here is 
valid for a 1 so that both expressions should be compared in the limit C^>1 where 
they indeed give the same result. Finally, when a < 0, the systems shrinks at long 
time and according to (36) the connective constant vanishes. 



The author enjoyed long discussions with Ingo Peschel and Ferenc Igloi about critical 
phenomena in parabolic geometries. 
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Appendix 

As suggested in the introduction the shape of the system may be considered as a 
perturbation to its critical behaviour in infinite geometry characterized by the scaling 
field \jC which, according to (1), may be either relevant, marginal or irrelevant 
depending on the sign of its scaling dimension 1 — za. This allows us to write down a 
scaling ansatz for the probability distribution 

1\ t b x ~ za 



p[^c) = b p [y¥— ) (A1) 

where the scaling dimensions are those of the unperturbed fixed point. 
With b= Ct a and z = 2 for the directed walk, (Al) translates into 

p (*-4) = cMc^^) (A2) 

so that for an irrelevant perturbation, a crossover towards the unperturbed critical 
behaviour occurs at t* =(7 2 /( 1_2a ) and a comparison with (7) shows that, in this case, 
the scaling function behaves asymptotically as a Gaussian 

/( "'" ) = 7S exp (-£> (A3) 

Furthermore (A2) can be inserted into (9) to get the asymptotic behaviour of the 
probability distribution in the case of a relevant perturbation. The scaling function 
then satisfies 

{1 -2a)v--af=-v^ + au- (A4) 

When a<l/2, v grows in time and the leading behaviour is obtained by keeping the 
two first terms on both sides to get 

which becomes exact in the strip geometry when a = 0. With f ~t(;(u)4>(v), ip even in 
u and vanishing for u = l, we find 

f(u, v)=J2 exp ( ~YZY ™<^mu) (A6a) 

m=0 ^ a ' 



A m =(m+0^= m = 0,l,2... (A6b) 

Together with (A2) this gives a properly normalized exact expression when a = 
whereas the term m = in the eigenvalue expansion provides the leading contribution 
(10) to the probability distribution. 
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